#ifndef RUNGEKUTTA_H
#define RUNGEKUTTA_H
#include "Integracao.h"

#ifdef RUNGEKUTTA_MINATO

class RungeKutta4ordem : public Integracao
{
    public:
            RungeKutta4ordem(){}
            RungeKutta4ordem(MODEL*m)
            { esferas=&m->SPHList; pModel=m;}
            void integra (SPHERE*);
            /*
            void integraPasso0 (SPHERE*);
            void integraPasso1 (SPHERE*);
            void integraPasso2 (SPHERE*);
            void integraPasso3 (SPHERE*);
            char step;*/
            void calcular();
            //double estimaVelocidade(double forca, double massa);
            //double estimaPosicao(double velocidade);
    protected:
            REAL H;
    private:

};
class RungeKutta4ordemParalelo : public RungeKutta4ordem
{
    public:
            RungeKutta4ordemParalelo(MODEL*m)
            {
                esferas=&m->SPHList; pModel=m; //step=0;
            }
            void calcular();
    protected:
    private:
};
#endif //RUNGE_MINATO


#endif // RUNGEKUTTA_H
